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Abstract 

Dynamic behavior of a weightless rod with a point mass sliding along the rod axis 
according to periodic law is studied. This is the simplest model of child's swing. 
Melnikov's analysis is carried out to find bifurcations of homoclinic, subharmonic 
oscillatory, and subharmonic rotational orbits. For the analysis of superharmonic 
rotational orbits the averaging method is used. The analytical results are compared 
with numerical simulation results. 



1 Introduction 



Oscillations of a pendulum with periodically varying length (PPVL) is one 
of the classical problems in mechanics, see [1-10]. We represent PPVL as 
a weightless rod with a point mass sliding along the rod axis according to 
periodic law. This is also a simple model of child's swing. In works [1,6-9] 
small oscillations of PPVL were studied. In [5,8,9] the instability domains of 
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the vertical position were found both analytically and numerically. Regular 
rotations and chaotic regimes were also investigated in [6,8,9]. 

In the literature oscillatory and rotational approximate solutions for PPVL 
were obtained with quasi-linear approach, where nonlinearity was assumed to 
be small as well as the excitation amplitude. The only exception is [10], where 
stable uniform rotations were found in the case of zero damping and special 
excitations. 

In the present paper we analyze boundaries in the parameter space for oscil- 
latory, rotational, and more complex (rotational-oscillatory, chaotic) regimes 
of PPVL assuming arbitrary (not small) nonlinearity. In Section 2 main equa- 
tions of motion of a pendulum with variable length are derived and given in 
non-dimensional form. In Section 3 we use Melnikov's analysis [11] to find 
bifurcations of homoclinic, oscillatory and subharmonic rotational orbits of 
PPVL. Melnikov's functions are obtained and compared with the results of 
numerical simulation. Similar analysis was done for a pendulum with vibrat- 
ing pivot in [14]. In Section 4 we find bifurcations of superharmonic rota- 
tional orbits with the use of the averaging method and compare them with 
the numerical simulation results. The methodological novelty in Section 4 is 
that the unperturbed system preserves the angular momentum (rather than 
Hamiltonian as in the Melnikov's analysis) and allows for non-small excitation 
amplitudes. 



2 Main relations 

Equation for motion of the PPVL is derived with the use of angular momentum 
alteration theorem and taking into account linear damping forces, see [5,6,8] 



where m is the mass, / is the length, 6 is the angle of the pendulum deviation 
from the vertical position, 7 is the damping coefficient, and g is the acceleration 
due to gravity. 

It is assumed that the length of the pendulum changes according to the peri- 
odic law 



where Iq is the mean pendulum length, a and Q are the amplitude and fre- 
quency of the excitation, <£>(t) is the smooth zero mean periodic function with 





I — l + aip(Qt), 



(2) 



2 



period 2ir. 



We introduce new time r 



Qt and three dimensionless parameters 




(3) 




where flo = ^Jf^ is the eigenfrequency of the pendulum with constant length 

I = Iq and zero damping. In this notations equation (1) takes the form 



((1 + eif(r)) 2 6)-+Pu(l + e V (r)) 2 6 + (1 + e<p(r)) u 2 sin(0) = 0, (4) 



where the upper dot denotes differentiation with respect to new time r. 



3 Melnikov's method: Perturbation of a Hamiltonian System 



Assuming 1 + £<f(r) > 0, equation (4) can be written in the following form 



Coefficients of nonlinear equation (5) explicitly depend on periodic function 
(/?(t) and three independent dimensionless parameters: the relative excitation 
amplitude e, the damping j3, and the inverse relative frequency of excitation 

OJ. 

Let us assume that parameters of excitation amplitude e and damping (3 are 
small of the same order, £ ~ /5< 1. Thus, we can say that dynamics of PPVL 
is described by the perturbed Hamiltonian system 




(5) 
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OH 

dv 



(6) 



v 



— + g 1 (8,v,r) + o(e), 



(7) 



where the perturbation function 



g±(9, v, r) = (2e sin(r) — (3cu) v + ecu 2 cos(r) sin(6 ) ) 
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is of the first order of smallness, i.e. gi(9, v, r) = 0(e). The following function 
H= ^--u 2 cos(9) (9) 



is the Hamiltonian of the unperturbed system 



9 = v, 

v = -co 2 sin(0), 



(10) 
(11) 



which is system (6)-(8) with e = and (3 = 0. The unperturbed system (10)- 
(11) describes the oscillations of the pendulum with constant length and zero 
damping, so that the integral H = const. 

Unperturbed system (6)-(8) has an oscillatory solution if H < uj 2 and a ro- 
tational solution if H > u 2 . If H = uj 2 the solution is a separatrix dividing 
oscillatory and rotational domains in phase space. 



3.1 Homoclinic bifurcations 



In order to apply Melnikov's criterion [11,12] to a homoclinic orbit we find the 
separatrix of the unperturbed system (10)-(11) that goes through the saddle- 
node point 9 — 7r, v — 0. For this point we have H = uj 2 and with the use 
of (9) obtain v 2 = 2u 2 (l + cos(#)) = 4u 2 cos(6 l /2). So, the separatrix has the 
following form 

v = 9 = ±2wcos y- J . (12) 



This equation allows for separation of variables and the following integration 
/ l + sin(g/2) \ 

111 1 cos(g/2) )=±*iT-'*)> ( 13 ) 



where tq is a constant of integration. Potentiation of (13) with some transfor- 
mations yields 

cos ( - ) = — — ; rr-, sin f - ) = tanh (to (r — r )) . (14) 

\2j cosh (u (r- r ))' \2j K J 



Melnikov's distance between stable and unstable perturbed separatrices is 
given by the following integral 
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oo 



M ± = J — gi (9,v,T)dr 

(2s sin(r) — j3u) v 2 (t) + su 2 v(t) cos(r) sin(6>(r))) dr 
which with the use of (12) and (14) takes the form 



-oc 



sin(r) dr ,„ , / dr 



M ± 8 2 / sm l r J dr 4/3 3 / 

J cosh 2 (a; (r — Tn)) J 



-oo 

oo 



cosh (u (r — Tq)) J cosh (lo (r — r )) 



' cos(r) t h( ^ (r ~: o)) dT, 

cosh (oo (r — r )) 



and has the following expression (see Appendix A) 



M± = 6 ! r £ sh 1 (ro)_ 
sinh 

This is a sign-changing function if the following inequality is satisfied 



e 4co 
— > 



sinhf —) . (16) 
{3 3tt \2u) k ' 

According to Melnikov's criterion [11] when condition (16) is satisfied then 
homoclinic structure and complex dynamics appear near the separatrix. Con- 
dition (16) means that complex dynamics appears only when the amplitude 
of excitation is sufficiently large with respect to the damping coefficient. The 
domain of possible complex dynamics defined by inequality (16) is depicted 
in Fig 1 and compared with numerical simulations in Fig 2 for /3 = 0.05. The 
minimum of the right-hand side of (16) is reached at u w 0.82 and can be 
evaluated as 0.948. Note that condition (16) is similar to that of oscillator 
with quadratic nonlinearity and external periodic excitation [13]. Similar in- 
equality for the pendulum with vertically vibrating pivot was obtained in [14], 
where the right-hand side function is three times greater than that in (16). 



3.2 Subharmonic bifurcations of oscillatory orbits 



In order to apply Melnikov's criterion to oscillatory orbits we find oscillatory 
solutions of the unperturbed system (10)-(11). For this reason we introduce the 
amplitude of oscillations A so that from the first integral for v = we have H = 



5 



— iJ 2 cos(A) and with the use of (9) we obtain v 2 = 2a; 2 (cos(0) — cos(A)) 
Auj 2 (sin 2 (A/2) - sin 2 (6>/2)). Thus, we can write the following equation 



9 = 2wsign(#) 



\ 



sin 



A 



— sin 



(17) 



which allows for separation of variables. In order to integrate (17) one usually 
introduces a monotonically increasing phase ip such that 




= k sin (ip) , 



where new variable k = sin (^J = y H ^j2 * s cane d modulus in elliptic inte- 
grals. Skipping intermediate transformations (see [15]) we have ip = am(c<;(r — 
To),k) and consequently 

v = 2uk cn (a; (r — r ) , k) (18) 



and 




dn(w (r - T ) , k) 




k sn(a; (r — r ) , k) 



(19) 



where am, dn, cn, and sn are the Jacobi functions. The k value follows from 
the resonance condition stating that period of oscillation and period of 
excitation 2ir should be in rational relation 

AK(k) , . 

^p = 2vrg, (20) 

u 



where p and q are relatively prime natural numbers and K(k) is the complete 
elliptic integral of the first kind. Thus, for oscillation motion in resonance 
p\q we have the following subharmonic Melnikov's distance with the use of 
expressions (18), (19) 



dH 



— 9i (9,v,t) dr 



o 

2-rrg 

(2esin(r) - (3uj)v 2 (t) + euj 2 v(t) cos(r) sin(0(r))) dr 
o 
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Fig. 1. Subharmonic bifurcation functions (e//5) which is the boundary in (23) for 
q = 2, q = 4, q = 6, and p = 1 converge to the homoclinic bifurcation function that 
approximates the domain of complex dynamics in (16). 

2nq 

= 4cu 2 k 2 J (2s sin(r) — (3ui) cn 2 (w (r — r ) , k) dr (21) 
o 

+ Aeu 3 k 2 J cos(r) sn(u (r — r ) , k) dn(w (r — r ) , k) cn(u; (r — r ) , k) dr. 
o 

Taking the integrals in Appendix B for p = 1 and even q — 2, 4, 6, ... we have 
the following Melnikov's distance 



MlM - ^ (^^Bmh ~ 40 ( E(k) ~ {k ' 2)im ^ ■ (22) 

where k' 2 = 1 — k 2 . This is a sign-changing function if the following inequality 
is satisfied 

£ > g ( BW _ (Otf(*W*MV (23) 



Condition (23) means that corresponding l:q resonant oscillations appear only 
when the amplitude of excitation is sufficiently large with respect to the damp- 
ing coefficient. The domain of possible oscillations denned by inequality (23) 
is depicted in Fig 1 and compared with numerical simulations in Fig 2 for 
(3 = 0.05. Note that condition (23) is similar to that of oscillator with quadratic 
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Fig. 2. Periods of oscillations divided by period of excitation (p : q) are shown 
with different colors on the parameter plane (lj, e) at (3 = 0.05. The correspondence 
between the colors and relative periods is shown by the color bar on the right, 
where stationary and rotational regimes are shown with white color, while chaotic 
regimes - with black. Subharmonic bifurcation functions for q = 2 (solid line) , q = 4 
(dot-and-dash line), q = 6 (dotted line), and p = 1 converge to the homoclinic bifur- 
cation function (bold solid line) and approximate the domains where corresponding 
resonant oscillations can occur, shown with blue (q = 2), purple (q = 4), and red 
(q = 6) colors. 

nonlinearity and external periodic excitation [13]. Similar inequality for the 
pendulum with vertically vibrating pivot was obtained in [14], where the right- 
hand side function is three times greater than that in (23). 



3. 3 Subharmonic bifurcations of rotational orbits 



In order to apply Melnikov's criterion to a rotational orbit we find the solution 
of the unperturbed system (10)— (11) for H > u 2 . Thus, with the use of (9) 



we obtain 9 = ±V2^jH + uj 2 cos(0) = ±2uk^l - sin 2 (#/2)/£; 2 , where "±" 
represents counter- and clockwise rotations. Since the model is symmetric with 
respect to the vertical axis, we will consider only the counterclockwise rotation 
("+" instead of "±"). This solution has the form 9 = 2am(cu/c (r — tq) , 1/k), 
so we have 



9 = 2cukdn[ujk(T - r ) ,-) , (24) 
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cos (^ J = cn(uk(r - t ) , ^ , sin^j = sn(uk (t - r ) . (25) 



The k value follows from the resonance condition stating that period of rota- 
tion 2K ( 1 J k '> and period of excitation 2% should be in rational relation 



ujk 



r = 2nq, (26) 



where r and q are relatively prime natural numbers. Thus, for rotational mo- 
tion in resonance r:q we have the following subharmonic Melnikov's distance 



OH 



M9/r = J ~fa 91{e ' v,T) dr 



2nq 

'(2e sin(r) - (3co) v 2 (t) + sco 2 v(t) cos(r) sin(0(r)) ) dr 
■ 8eco 2 k 2 J (2e sin(r) - dn 2 (uk (r - r ) , dr 



o 

2vn} 





+ 4eo; 3 /c ^ cos(r) sn^uk (r — r ) , — ^ dn^wA; (r — r ) , — ^ cn^wA; (r — To) , — ) 
o 



Taking the integrals in Appendix C for r = 1 and q — 1, 2, 3, ... we have the 
following Melnikov's distance 



^-^fesgr*"®). (27) 



where IT = K(Jl-l/k 2 ). This is a sign-changing function if the following 
inequality is satisfied 



e Au 2 k 

p > i^r \k 



Condition (28) means that corresponding l:q resonant rotations appear only 
when the amplitude of excitation is sufficiently large with respect to the damp- 
ing coefficient. The domain of possible rotations defined by inequality (28) is 
depicted in Fig. 3 and compared with numerical simulations in Fig. 4 for 
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Fig. 3. Subharmonic bifurcation functions (e//3) which is the boundary in (28) for 
q = 1, q = 2, q = 3, and r = 1 converge to the homoclinic bifurcation function that 
approximates the domain of complex dynamics in (16). 

= 0.05. Although numerically found rotations 1:2 are not monotone. These 
are rotation-oscillation regimes. Note that condition (28) is similar to that 
of oscillator with quadratic nonlinearity and external periodic excitation [13]. 
Similar inequality for the pendulum with vertically vibrating pivot was ob- 
tained in [14], where the right-hand side function is three times greater than 
that in (28). 



4 Averaging method: Superharmonic bifurcations of rotational or- 
bits 



We study resonant rotations r:q = 1:1,2:1,3:1, . . . under assumption that uj 
is small (e.g., because of the small gravitation g or high excitation frequency 
Q) of the same order as damping, ~ oj <C 1, while the excitation amplitude 
e is not small. Then, unperturbed system (equation (4) with (3 = u = 0) 
has the specific angular momentum or sector velocity s = (1 + e(p(r)) 2 6 as its 
first integral. The unperturbed system is not Hamiltonian, so we cannot apply 
Melnikov's analysis to obtain domains of existence of corresponding rotations. 
For this purpose we use the method of averaging, [2,17]. Equation (4) can be 
written in the form of the system 
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homoclinic 
r:q=1:1 rotations 
r:q=1 :2 rotations 
r:q=1 :3 rotations 




1:2 



0.5 0.6 0.7 O.f 
relative frequency <o 

Fig. 4. Periods of regular rotations divided by period of excitation (this ratio is r : q 
for oscillations) are shown with different colors on the parameter plane (w, e) at 
/3 = 0.05. The correspondence between the colors and relative periods q is shown 
by the color bar on the right, where stationary and chaotic regimes are shown with 
white color. Subharmonic bifurcation functions for q = 1, q = 2, q = 3, and r = 1 
converge to the homoclinic bifurcation function and approximate the domains where 
corresponding resonant rotations can occur. 
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a+^(r)) 2 ' 

s = f(0, s,t), 



(29) 
(30) 



where the perturbation function 



f(6, s, t) = -f3us - (1 + £<p(r)) uj 2 sin(0) 



(31) 



is of the second order of smallness, meaning f(6,p,r) = 0(u 2 ). The unper- 
turbed system has the solution 



dr] 



(l + etp(ri)y 



(32) 



where $(r, e) = Jq ^ 1+ f^ v ^2 and is the constant phase shift. In the first 
approximation we can derive the averaged equation of (30) substituting in it 



11 



expression (32) and taking the average of its right-hand side, see [17] 

2nqr 

2nqr 
* o 



£=-!-/ f(s^(r,e)+^s,r)dr (33) 



where the integral is taken as if corresponding averaged variables s and $ 
are constant. The steady state value s follows from the resonance condition 
written with the use of solution (32) and stating that period of rotation and 
period of excitation should be in rational relation 



2wq 

dr 



^r = s I — = sg$(2vr,e), (34) 



o 



:i+^(r)r 



where r and q are relatively prime natural numbers. Steady state (s = 0) value 
of $ we find from (33). 1 Thus, with the use of (34), for rotational motion in 
resonance r:q we have 



(27rr) 2 /3 r r $(r,e) 



/ (1 + eip(r)) sin ( 2ir - x 
W $(2tt,£) J 1 ^ V JJ \ g$(27r,e) 

= A(e, r, q) cos^J + B(e, r, q) sin^J (35) 
where A(e,r,q) and B(e,r,q) denote the following integrals 



+ dr 



^, r , ?)= /( 1+ e^))sin(2.:-ifc|)dx, 

2nrq . t \ \ 

^^)=/(l + s^))cos(2^^) dr. 



(36) 
(37) 



1 One can also write the differential equation for slow phase shift as •& = 
<&(t, e) /(s <E>(t, e) + s,r), but it has unbounded right hand side because of the 
growing multiplier $(r, e) ~ r. That is why this equation does not fit for standard 
Bogolubov's method of averaging, see [2]. Thus, we can obtain here only the first or- 
der approximation with the use of averaged equation (33) . In order to obtain higher 
order approximations the general averaging scheme by Volosov can be applied to 
system (29)-(30), see [17]. 
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Fig. 5. Superharmonic bifurcation functions (to/ (3) which is the boundary in (38) 
for r = 1, r = 2, r = 3, and q = 1. 

Thus, the domain of existence of corresponding regular rotations is approxi- 
mated by the following condition that equation (35) has a solution 



which right hand side can be calculated numerically for any particular e, r, 
and q, see Fig. 5, where ip(r) = cos(r). 

In order to find the boundaries for rotational regimes r = 1, 2, 3 the first ap- 
proximation is enough because excitation amplitude e is not small in contrast 
to the quasilinear approach in [9], where both e and u are assumed to be 
small and higher order approximations of the averaging method are needed to 
obtain similar boundaries. 



5 Conclusion 

For the pendulum with variable length we derived analytical formulas for the 
boundaries of bifurcations in the space of three parameters: relative frequency 
of excitation, amplitude of excitation, and damping. The boundaries for ho- 
moclinic bifurcation separating the domain of only stationary and oscillatory 





P $(27r,e) v /A 2 (£,r,g) + J B 2 (e,r, qj' 



13 




relative frequency <o 

Fig. 6. Periods of regular rotations divided by period of excitation (this ratio is r : q 
for oscillations) are shown with different colors on the parameter plane (w, e) at 
/3 = 0.05. The correspondence between the colors and relative periods r is shown 
by the color bar on the right, where stationary, oscillatory, and chaotic regimes 
are shown with white color. Superharmonic bifurcation functions for r = 1, r = 2, 
r = 3, and r = 1 bound from below corresponding regimes obtained numerically. 

regimes from the domain of more complex dynamics, subharmonic oscilla- 
tions and subharmonic rotations are obtained using Melnikov's method under 
assumption of small damping and excitation amplitude. For the analysis of su- 
perharmonic bifurcations of rotational orbits the method of averaging is used 
assuming smallness of relative excitation frequency rather than that of excita- 
tion amplitude. Both methods allow to obtain in the first approximation the 
basic rotational orbit with angular velocity equal to the excitation frequency, 
1:1. Small excitation frequency (or small gravity) allows to introduce the un- 
perturbed system with the conservation of angular momentum, so that faster 
rotations are found in the first approximation by the method of averaging. 
In Figs. 2, 4, and 6 it is shown that the boundaries for complex dynamics, 
subharmonic oscillations, and superharmonic rotations are in good agreement 
with the results of numerical simulation. 
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Appendices 



A Melnikov function for homoclinic orbit 



M d 



8eo; 2 Ji — A(3u 3 l2 + Aeu 3 ^, where the first integral 



sin(r) dr 



cosh (u 



T 



TO)) 



1 r sin(r + tj/uj) dr] 
uj J 



cosh (77) 



sin(r ) f cos(r]/uj) dr] cos(r ) f sm(r)/u) dr] 



cosh 2 (77) 



uj 



cosh 2 (77) 



(A.l) 



The integral S1 ^^^ r? is zero because its integrand is an odd function, 
while the other integral has even integrand and can be calculated as follows 
•/-00 C °osh^ ^ = u sinh(7r/2^) ' nence 5 the first term has the expression 



h 



7rsin(r ) 



uj 2 sinh(7r/2a;) 
The second integral can be calculated as follows 



(A.2) 



dr 



cosh (a; (r — r )) 



ds 



cosh (s) cu 



(A.3) 



while the third integral can be converted to the first one via integration by 
parts using the relation 



sinh(s) ds 
coshes) 



-id 



1 



2 cosh (s) 



h = 



UJ 



cos(r) sinh {u (r — r )) 
cosh 3 (a; (r — r )) 



cos(r + tj/uj) sinh (77) 
cosh 3 (77) 



dr, 



dr] 



1 cos(r + r]/ui) 



2uj 2 cosh (r/) 



2uj< 



sin(r + rj/u) 
cosh 2 (77) 



*, = -£. (A.4) 
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Thus we have Melnikov function M ± = 8euj 2 h-A^uj 3 U+Aeuj 3 I 3 = 6w lf n t\ 

1 " z sinh(7r/2ai) 

8f3u 2 . 



B Melnikov function for subharmonic oscillations 



M«/p = 8eu 2 k 2 h - A/3u 3 k 2 I 2 + Aeu 3 k 2 I 3 , where 



2-rrq 



I\— J sin(r) cn 2 (a; (r — r ) , k) dr 





2irq 

2wJ 3 , 
d ' 



= — cos(r) cn 2 (w (r — r ) , fc) 
= -2coI 3 , (B.l) 

7 2 =i cn 2 (a;(r-ro),A;) dr, (B.2) 


2irq 

7 3 = J cos(r) cn(u; (r — r ) , fc) sn(u; (r — r ) , fc) dn(a; (r — r ) , k) dfp.3) 


where we use the formula dc ^f > = — sn(u) dn(u). Thus, we have M q / p = 
-12eu 3 k 2 I 3 - A/3u 3 k 2 I 2 . 

2rcq 

h— J sin(r) cn 2 (u; (r — r ) , k) dr 


1 r 

= — / sin(r + s/oj) cn 2 (s, k) ds 



= snHroj /■ cos ( s y cn 2( S) n ds / b _ 4 n 



The integral J3 vanishes except for p = 1 and even g. In this case we have 
from (312.02) in [16] 

h = % (m - (k ,2 )K(k) 

71 sin To 



k 2 u smh(K(k')) ' 

where k' 2 = 1-k 2 . So we have M"/ 1 = 4w 4 (gff sin ^ fc/)) - 4/3 (£ , (/c) - (k' 2 )K(k 
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C Melnikov function for subharmonic rotations 



Mi/ r = %eu 2 k 2 h - Af3u 3 k 2 I 2 + 4eu 3 kl 3 



2nq 

I\ — J sin(r) dn 2 (^uk (r — r ) , — \ dr 



dn 2 ( uk (r — r ) , — 



2 " 2uj i - 2uj i (C 1) 

o "T /3 -"T /3 ' (CU) 



I 2 = J dn 2 (U;(T-7o),£) dr, (C.2) 
o 

2nq 111 

h = J cos(r) sn^wA; (r - r ) , -J dn^wA; (r - r ) , cn\uik (r - r ) , ^(JCdS") 



where we use the formula ddn (™M k ) = — sn(w, 1/k) cn(w, l//c)/fc 2 . Thus, we 
have Mil r = -12eoo 3 k 2 I 3 - 4(3(Jk 2 I 2 . 

The integral /3 vanishes except for r = 1. In this case we have 

2 uk (&) ' 
7T sin Tq 



k 2 u sinh(i^') 



where K' = K[ x /l - 1/k 2 ). So we have M q/1 = -2uj 2 ( ^^fr - W(| 
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